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Abstract 

Using the MIT bag model, we calculate the free energy of droplets of quark- 
gluon plasma in a bulk hadronic medium, and of hadronic bubbles in a bulk 
quark-gluon plasma, under the assumption of vanishing chemical potentials. 
We investigate the validity of the multiple reflection expansion approxima- 
tion, and we advise a novel procedure for calculating finite-size corrections to 
the free energy of hadronic bubbles in a bulk quark-gluon plasma. While our 
results agree largely with earlier calculations, we show that the usual multiple 
reflection expansion should be used with caution, and we propose a modifi- 
cation of the multiple reflection expansion, which makes this approximation 
agree nicely with direct numerical calculations. The results should be of rele- 
vance in connection with the cosmological quark-hadron transition as well as 
for ultrarelativistic heavy ion collisions. 
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I. INTRODUCTION 



The quark-hadron phase transition is of significant interest in connection with ultrarel- 
ativistic heavy ion collision experiments, the interior of neutron stars, and the evolution of 
the early Universe. A calculation from first principles using QCD is at present impossible, 
but lattice-QCD studies have shed some light on the transition, for instance demonstrating, 
that the transition is apparently first order for pure glue, whereas the order for physical 
QCD is still a matter of investigation. 

Awaiting more definite answers to come from such investigations, numerous studies have 
been performed using phenomenological models, in order to gain insight into the physics 
of the transition. Many such studies have used the MIT bag model, which in a relatively 
simple manner incorporates confinement in terms of a set of boundary conditions for quarks 
and gluons. 

A very interesting result of a detailed study within the MIT bag model was presented 
by Mardor and Svetitsky jl], who considered the zero chemical potential case of relevance 
for the cosmological quark-hadron transition. For a droplet of quark-gluon plasma within 
a bulk medium of pions, a direct numerical calculation of the partition sum using quark 
and gluon energy levels led to a behavior of free energy as a function of radius, F(R), as 
expected for a first order transition, namely a minimum of F for R = when T is below 
the transition temperature, T , and an energy barrier for R of order a few fm separating a 
local minimum at R = from the true minimum (diverging negative energy) for R — > oo. 

To treat the "inverse" problem of a vacuum (hadron) bubble within a bulk phase of 
quark-gluon plasma, the authors employed a phase shift formula to calculate the changes in 
quark and gluon density of states stemming from the presence of the hadron bubble; again 
calculating the contribution to the free energy by a direct numerical integration. In this 
case, a peculiar feature was observed, namely that F(R) had a negative minimum for radii 
of 1-2 fm, even for T > To, apparently indicating an instability of the quark-gluon plasma 
above To, since there was no energy barrier to prevent formation of hadron bubbles. 

An interpretation of the result was put forward in terms of an expansion of the free 
energy in terms of volume, surface, and curvature contributions, 

4 

F(R) = AP-ttR 3 + a47rR 2 + a8irR + ... (1) 
3 

Here R is the radius of the droplet or bubble, AP is the pressure difference between quark 
and hadron phases (with AP = defining the transition temperature, To), a is the surface 
tension, and a the curvature coefficient, where volume, surface and curvature terms can be 
calculated from the smoothed quark and gluon densities of state within the MIT bag model 
(see below). The results of Ref. |1J were apparently well reproduced under the assumption 
that a vacuum bubble behaves like a plasma droplet turned inside out, so that the radius 
changes sign. This leaves the area term unchanged, but volume and curvature terms change 
sign. Or, if R is defined to be positive in the expression above, then a is unchanged, but 
pressure difference as well as curvature coefficient changes sign when going from the case of a 
plasma droplet to that of a vacuum bubble. The fact that the curvature contributions from 
massless quarks and gluons are very large compared to the surface contributions coming 
only from the massive s-quarks, could explain that an energy barrier for the plasma droplet 
could turn into an energy minimum in the reverse case of a hadron bubble. 
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The authors of Ref . |1] were careful to point out a number of reasons to be cautious about 
the result. First of all that the MIT bag model is clearly just a phenomenological model, and 
also that the radii of relevance for the interesting bubbles/droplets were perhaps too small 
to justify the ideal gas approximations. But if the result was of a physical nature, it did have 
important implications for the understanding of the quark-hadron transition . And the 
procedure of R — > —R gave a simple recipe for treating other situations, such as quark- 
hadron mixed phases in neutron stars. Some consequences, however, were rather strange. 
For instance, it apparently pays energetically to fill a strangelet with vacuum bubbles, so 
that it looks more like a Swiss cheese than like a uniform mixture of quarks []5| . 

The aim of the present paper is to compare several different calculations of the free energy 
of a vacuum bubble embedded in quark-gluon plasma as well as a quark-gluon plasma droplet 
within a bulk phase of hadrons. For the plasma droplet we focus on a direct sum over states 
compared with a multiple reflection expansion, showing that terms beyond volume, surface, 
and curvature in the free energy are necessary to avoid unphysical behavior for small radii. 
We demonstrate how the next important contributions to F (proportional to T\n(RT), and 
to T) arise naturally if the density of states is truncated below some value of kR (where k 
denotes momentum) instead of integrating over unphysical, negative values for the density 
of states all the way from k = 0. For the vacuum bubble we confirm the results of the 
phase shift approach of Ref. |l| by comparing to a more direct sum over states approach 
introduced below, which we refer to as the concentric spheres method. Again we show how 
an improvement of the multiple reflection expansion leads to correction terms in F, such 
terms arising naturally from a truncation of the density of states. 

The general framework and basic equations are described in Section [H]. In Section IT 1 1 
we present our numerical results, largely confirming the calculations of Ref. Our results 
show how and why the usual version of the multiple reflection expansion is not always 
accurate. In Section |IV| we show how further terms in the analytical expansion of the free 
energy proposed in the literature improves the agreement with the numerical results, and 
we show how a physically motivated truncation of the density of states from the multiple 
reflection expansion resolves most of the problems encountered in Section III. Section 
contains our conclusions. 



II. THEORETICAL FRAMEWORK 

In this section, we give the basic equations needed for an analysis of the quark-hadron 
phase transition within the MIT bag model. We consider the case of zero chemical potential 
which is of particular relevance to the cosmological quark-hadron transition, but also of 
interest for ultrarelativistic heavy ion collisions. The quark-gluon plasma is taken to consist 
of three quark flavors (u,d and s), the corresponding antiquarks, and eight non- interacting 
gluons, these particles being described by the MIT bag model presented below. The hadron 
phase is considered a mixture of the 3 pions n ,^, since all other (much heavier) hadrons 
contribute only insignificantly to the free energy. Further, we shall assume that the pions 
contribute volume terms only (see Sec. [II B|) . We have taken m u = ma = 0, m s = 150 MeV 
and m n = 138 MeV. 
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A. The MIT bag model 



The MIT bag model ||[7] is defined by the Lagrangian 

L= f d 3 x(C QCD -B). (2) 

£-qcd is the usual QCD Lagrangian density, and L = outside the bag. Q is the bag-volume. 
B > is a phenomenological parameter, the bag constant, which models the difference in 
energy density between the perturbative vacuum inside the bag and the non-perturbative 
QCD vacuum outside the bag. Requiring the action W = J* t 2 dt L to be stationary with 
respect to variations of the fields yields the equations of motion. 

At the surface of the bag, the fields are taken to satisfy boundary conditions which 
correspond to the fields being confined inside the bag-volume. We neglect gluon-exchange 
interactions. 

The equations of motion for the fields become the Dirac equation for the quark fields, 
and the source-free Maxwell equations for the gluon fields. The complete set of equations 
governing the behavior of the fields, including the boundary conditions, is 

(17"^ - m)V(x) = 0, xett (3) 
d li F'"(x) = Q, xeQ (4) 
in l sy IJ "$(x) = x e dVt (5) 

n M F^(x) = 0, xedVl (6) 

in the following notation: x = (x°,x) is a space-time four-vector, dQ is the surface of 
the bag-volume Q, for x G dfl we define n^(x) = (0,—x/\x\) as an inward-directed unit- 
normal three- vector to the surface of the bag, *$?(x) is the quark-spinor (there will be one 
for each quarkflavor (u,d,s,...) and one for each of the three color states of a quark) and 
F^ u (x) = d^A u (x) — d u A fl (x) is the (non-interacting) gluon field (there are eight copies of 
this field). 

We fix the bag constant by demanding bulk pressure balance at the transition tempera- 
ture. Somewhat symbolically the bag constant is thus determined by the equation 

73 j. J ^Fquarks 9Fgi uons 9Fpi ons 1 , . 

~ v^L \ dv dv dv j T=To ' 1 ' 

F being the free energy. In the following, we shall set the transition temperature to T = 
150 MeV, thus fixing the bag constant B = 312.6 MeV/fm 3 = (221.4 MeV) 4 . 

We can immediately write down the expression for the gluon field, since this is just 
the solution to the source-free Maxwell equations. Expressing the gluon field in terms of 
color-electric and color-magnetic fields, writing A^{x) = (V(x),A(x)), V x A = B(x) and 

-Vy(i) — 9 q^o > = E(x), there are two sets of solutions to (|]), labeled TM and TE (the 
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I = fields are absent since for 1 = 0, the only solution to the source- free Maxwell equations 
is Boo = E m = !): 

K^-^«V-£^;; w (8) 



where 



Bfr? = -i fi(kr) x x (VY lm {e, 0)) (10) 
= -iVx(f/,(Ax)Y, m (M)), (11) 



e4 vx tt) (12) 

= ly x V x (f /,(At) Y lm (0, </>)), (13) 



C = ~i fi(kr) x x (Vy, TO (0, 0)) (14) 
= -iVx(f/,(A;r)Y, ro (0,0)), (15) 



C £ = -JVx&) (16) 
= ~ V x V x (f /,(A;r) y, m (0, <j>)), (17) 

Yi m (6,(j)) are the usual spherical harmonics, /j(^) = a>ji{z) +bni(z), the spherical Bessel- 
functions ji(z) and 71/(2;) being the two linearly independent solutions of the equation 

z 2 g"(z) + 2zg' (z) + {z 2 - 1(1 + l))g(z) = 0, 

r = \x\, and = (\k\,k). The constants a and b appearing in the function fi, and the 
possible values of k, must be fixed from the boundary conditions (|6|). Expressed in terms of 
the fields E and B, these boundary conditions read: 

x-E = xxB = 0, xedQ. (18) 

a. Extension of the MIT bag model. The MIT bag is a finite region of space(-time) to 
which quarks and gluons are confined by boundary conditions (^)-(^) corresponding to no 
flux of plasma out of the droplet. We shall refer to this configuration as a plasma droplet. 
However, in the following we shall also use the MIT bag model in a slightly different way, 
namely the case where quarks and gluons are kept outside a finite region of space, cf. Fig. |1|. 
The equations describing the quarks and gluons in the second configuration, the "vacuum 
bubble", are still (|3|)-(|6|), but now using ri\x) = (0,x/\x\) in the boundary conditions. 
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B. Thermodynamical relations 



For a system of non-interacting fermions (upper sign) or bosons (lower sign) we can 
calculate the free energy for each particle degree of freedom as 

oo 

F(T, V) = tTJ2 ln (! ± e- E ' {v)/T ) (19) 
i=i 

where Ei(V) = \Jm? + k 2 {V). In the continuum case we have 

F(T,V) = tT J d 3 k p(k,V)\n(l ± e- Vm2+k2/T ) } (20) 

p(k, V) being the density of states, defined such that p(k, V)d 3 k is the number of states in 
the volume V with momentum in d 3 k around k. The importance of the free energy stems 
from the fact that the configuration realized in nature is characterized by a minimum in this 
free energy. 

In the following we shall speak of the volume part resp. the surface part of the free 
energy. In the case of non-interacting Dirac particles and non-interacting gluons (these are 
just Maxwell fields) which are the particle species relevant to us, the density of states p in 
any sufficiently large volume V contains a term proportional to the volume. In fact[], 

p(V^oo)~ 9i ^ (21) 

independent of which particle species we consider. Therefore, also the free energy will contain 
a term proportional to the volume of the system. We name this term the volume free energy. 
The total free energy being F tot , we can write F tot = f voi V + F sur , where F sur /V — > as 
V — >oo, and f vo i does not depend on the volume. We shall call F sur the surface part of the 
free energy, or simply the surface free energy. 



C. The multiple reflection expansion, MRE 

The multiple reflection expansion (MRE) is an approximation for the density of states, 
also commonly referred to as the asymptotic expansion of the density of states. Since we 
only consider systems with spherical symmetry, we define the spherically symmetric density 
of states p(k,V) = 4irk 2 p(k,V). Consider a spherical volume V = ^-R 3 of quarks and 
gluons, described by the bag model. The MRE for this system, as a sum of volume, area, 
and curvature contributions, valid for sufficiently large volumes, is 

Vk 2 

Pi(k,V) = — Y + f A , i (k/m)kA7iR 2 + fcAk/m)87iR+..., i = q,g (22) 

where 



1 gi accounts for spin (helicity) degeneracy. 
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1 2 

fA, q (k/m) = - — (1 arctan(fc/m)) (23) 

07T 7T 



1 3k ix 

fcMm) = ^(1 - — (- - arctan(/c/m))) (24) 



= (25) 



Here, A stands for area, C for curvature, and indices q and g denote quarks and gluons. 
Note that lim m ^ fc, q (k/m) = -^2, and that \im. m ^ f A , q {k/m) = 0. 

The MRE was developed by Balian and Bloch ||, and the above expressions for the 
area- and curvature-terms have appeared in the literature. The area term for quarks is 
given (though not derived) in ||10|| . The curvature term for massless quarks seems to appear 



explicitly for the first time in pT| , whereas the full expression (|4]) for massive quarks is 
introduced in Hl2]| . The gluon expressions, valid for non-interacting gluons, is calculated in 

As indicated by the dots in (|22"1), the expression for pi(k, V) should in principle contain 
terms proportional to 1/-R, 1/ R 2 etc., but as these terms become small in the limit of large 
R, and since the MRE is an approximation valid for large systems, these terms are usually 
neglected. However, we shall see in the following that the MRE as it stands in (^) is 
not only inaccurate, but also unphysical at small radii, having negative density of states. 
Further, we shall argue that, when used in calculations of the free energy, the MRE ( |2"2"D 
containing only area- and curvature-terms leads to errors even at larger radii, where the 
MRE itself is a good approximation to the density of states. We also suggest a solution to 
these problems. 

Everywhere in the following, unless explicitly stated, reference to the MRE means the 



approximation (|22|) to the density of states without further correction terms like 1/R, 1/R . 



D. The inverse multiple reflection expansion, MRE(-fi) 

Because of the "symmetry" between the two situations (i) quark-gluon plasma confined 
by MIT bag boundary conditions within a sphere of radius R (a "plasma droplet"), and (ii) 
quark-gluon plasma kept outside a sphere of radius R by MIT bag boundary conditions (a 
"vacuum bubble"), it has been argued |IJ that there should exist a simple relation between 
the density of states in the two cases, i.e. that the density of states of quarks and gluons 
in the case of a vacuum bubble can be found from the expressions fl2~2"|)-(p6|) for a plasma 
bubble by simply inverting the sign of R. We shall refer to this hypothesis as the MRE(— R). 
Since the MRE(— R) is derived from the MRE, we expect the MRE(— R) to have problems 
related to those of the MRE mentioned in the previous subsection. 
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E. The phase shift approach 



In this subsection we briefly describe the phase shift approach to calculating the free 
energy of a vacuum bubble. The phase shift formula ( p7f) was introduced in this context 
by Mardor and Svetitsky The phase shift approach is based on a relation between the 
density of states and the scattering phase shifts^: 

. .,. 1 d5i(k) .„„. 



For a derivation of this relation in the non-relativistic case, see e.g. Here, A/?; is the 
change in the density of states (at a given angular momentum) induced by the scatterer. 

In order to use this phase shift approach to calculate the free energy, we need the scat- 
tering phase shifts for quarks and gluons. These are derived from the defining equations ([|) 
and (H) for the quarks, and from fl8|), @ and ( fL8|) for the gluons. The phase shift for the 
j-component of the quark field (j = 1/2,3/2, ... is total angular momentum) is the sum of 
two components 

5j(k) = 6 l = J ~ 1/2 (k) + 6 l = j+1/2 (k), (28) 
where, for a surface of radius R, 

1 {k) - mtm {n,(kR)± 1 ^-n l UkR) ' <29) 



and E = yrn^ + k 2 , q = u, d, s. The phase shift for the gluon field also consists of two parts, 
Sf M (k) and 5f E (k), where 



5j E {k) = arctan ( ) (30) 



arctan 



' 3l (kR)(l + l)-kR 3l+l (kR) \ 

ni{kR){l + l)-kRni +1 {kR)) ' 1 ' 



and 



- — (SI) • < 32 » 

again for a surface of radius R. I = 1,2, ... labels orbital angular momentum, and here, as 
opposed to the quark situation, it is a good quantum number. 

Knowing the phase shifts, the contribution to the free energy from the quarks and gluons 
outside a vacuum bubble of radius R is calculated using (|20|), so that 



2 In the case of spherical symmetry, the phase shifts 5i(k) are defined such that the effect of the 
scatterer is to change the spatial part of the wave function far away from the scatterer for a given 
angular momentum I from oc -r- sin(/cr — Itt/2) to oc -r- sin(/cr — Itt/2 + 5i(k)). 
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Fi(T, R) = T9i- r dk ln(l ± e ~ W) (33) 

The label i stands for different particle types (quarks, gluons) and angular momentum. 
Again, the upper sign applies to the fermions (quarks), lower sign to bosons (gluons). The 
appropriate degeneracy factors are g qua rk = 6 and g g i ucm = 8. 

We make two remarks about the formulae fl2"T|) and ([33]). (i) The free energy ([33]) includes 
the contribution from the excluded volume, (ii) When the phase shifts contain functions with 
multiple branches, like the arctan function in our case, we choose the branch which makes 
the phase shifts continuous functions of the energy. (In the case of potential scattering where 
the potential obeys certain integrability conditions, one can prove that the phase shifts are 
continuous functions of the energy fl5fl .) 

1. The free energy in the limit RT — > 
By expanding the Bessel- functions appearing in (|29|), (|3T| ) and ( p2|) as 



oo 



Ji(x) = J2*k(l)x l+2k (34) 

fc=0 

and 

oo 

n l (x) = J2Ul)x 2k - l -\ (35) 

fc=0 

(valid for / > 0) and keeping only the lowest order terms, we obtain via (|33|) the follow- 
ing analytical expressions valid for RT 1 for the surface free energy of (one flavor of) 
massless quarks (index j and / means that we consider each angular momentum component 
separately) 



FURT) io (2j + l)(2^+ 2 -l) 



12 ^ ' ■ - V?F a q (j) (RT)^ +2 B J+3/2 (36) 



2j+2 



T (2j + 3) 

and for the eight gluons 

F s, g (RT) (2/ + l)(27rf+ 1 
T 21 + 2 



a 9 (l)(RT) 2l+l B l+1 (37) 



where 



r .x _ aoU + 1/2) bp(j + 1/2) - a (j - 1/2) bp(j - 1/2) 
aq[3) ~ b 2 {j + 1/2) ' [6 * } 



m _ ^0(0 flo(Q (j + 1) ,o Q X 

aa[l) - b (l) + b (l)(l + l)-b (l + l) [6J) 



and 
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(-l) fc 

ttk{l) = 2 k k\l-3-5---(2l + 2k+iy (40) 

h(l) = (-l) fe+1 2 *k\(l-2l)(3-2l)---(2k-l-2l)- (41) 

(The factors (1 — 2Z) (3 — 21)... in the denominator of fl4"T|) appear only when 2A; — 1 > 1). 
The _B n appearing in (j36|) and (|37| ) are the Bernoulli numbers, defined by 

2 4 P 



- 2 B '¥- ft ¥ + B 'ffi-- < 42) 
the first few of these being 

B x = 1/6, B 2 = 1/30, B 3 = 1/42, E 4 = 1/30. (43) 



On the basis of (|36|) and (37), we conclude that the first energy term of importance for 
R — > is proportional to R 3 ; no terms proportional to R or R 2 appear in this limit. This is 
in contrast to the MRE(— R) conjecture, where a curvature term proportional to R dominates 
for R — > 0. The difference is clearly demonstrated in the figures in the next section as a 
difference between zero and finite slope of the free energy for R — > 0. 



F. The concentric spheres method 

The concentric spheres method is a new way to calculate the surface contribution to 
the free energy of quarks and gluons outside a vacuum bubble. The idea is to extract 
the contribution from the inner surface to the total free energy of the concentric spheres 
configuration in Fig. £| Assuming that the splitting of the free energy in volume and surface 
contributions is valid, we can extract the free energy contribution from the inner surface 
(cf. Fig. |^) from a calculation of the total free energy of particles contained between two 
concentric spheres as 

F S urface{Rl) = Ftotal (Rl , Rl) + {Rx)-F total (R 2 ), (44) 

where F tota i(Rx, R 2 ) is the total free energy (including both surface contributions) of the 
particles contained between the spheres with radii R\ and R 2 , F vo i ume (Rx) is the volume 
contribution to the free energy of the particles in a sphere of radius Rx, and F tota i(R 2 ) is 
the total free energy (including surface contribution) of the particles in a sphere of radius 
R 2 . It is precisely F sur f ace (Rx), the surface free energy of the particles outside a sphere, that 
we are interested in. Here, the terms F tata i(Rx, R 2 ) and F total (R 2 ) are calculated directly 
by summation over energy levels, whereas F vo i ume (Rx) is calculated from the MRE (using a 
positive radius). Explicitly, the term F tota i(Rx, R 2 ) in the case of quarks of mass m q is 

00 

F quarks (Rx,R 2 ) = -6T ]T (2j + l) ]T JXl + e-^"^), Ejln = y/k] ln + m*, 

j=l/2,3/2 l=j±l/2 n 

(45) 
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where fcj,z=j±i/2,n is the n'th solution of the equation 

(aiQj^kRi) - ji{kR!)) (a{k)n l _ 1 {kR 2 )+n l {kR 2 )) - 
(a(k)ji-i(kR 2 ) +ji(kR 2 )) (a(k)n l ^ 1 (kR 1 ) - n l (kR 1 )) = (46) 
for I — j + 1/2, and 

(a(k)ni + i(kRi) + n^kR^) (ji(kR 2 ) - a(k)j l+1 (kR 2 )) + 
(a(k)n l+1 (kR 2 ) - ni {kR 2 )) (fcfli) + 3i(kRi)) = (47) 

for I — j — 1/2, with a(fc) = , fc 2 . In the case of gluons, the term F to tai(Ri, R2) is 

•w k -\-rriq -\-rriq 

00 

F gluons (R 1 ,R 2 ) = 8T]T(2Z + 1) E ln U " ^ KJT ), (48) 

Z=l a=TM,TE n 

where kf n is the 77,'th solution of the equation 

Jl (kR 2 ) mikRr) - jtikRx) rn{kR 2 ) = (49) 

for the TM gluons, and 

(/ + 1) 2 {jiikRt) m{kR 2 ) - j t {kR 2 ) mikRt)} + 

(I + 1) kR 1 {ji(kR 2 ) n l+l {kR{) - j l+1 {kR x ) ni (kR 2 )} + 
(I + 1) kR 2 {j l+1 (kR 2 ) m{kR x ) - UkR,) n l+l {kR 2 )} + 
kR\ kR 2 {ji +1 (kR 1 ) n l+l {kR 2 ) - ji +1 (kR 2 ) n l+1 (kR 1 )} = 0, (50) 

for the TE gluons. The term F tota i(R 2 ) is just the free energy of a quark- or gluon-droplet, 
so this is calculated using Eqs. (|53|)-(|57]). Finally, the volume term F vo i ume (Ri) is 

poo Vfc^ 

F volume {R l ) = T9iT dk-—ln(l±e~ E ^ T ),i = q,g (51) 
Jo 2ir z 

where V = ^f-R\, g q = 12 for each flavor, and g g = 16. 

Formally, the splitting in surface and volume terms is appropriate only when R 2 — R± — > 
00. However, calculations for 10 fm< R 2 < 20 fm suggest that the concentric spheres method 
yields the correct free energy contribution from the inner surface up to R± ~ R 2 /2. 



III. NUMERICAL RESULTS 



We are now going to use the different techniques described in the previous section to 
calculate the free energy of (i) a plasma droplet in a bulk hadronic medium, and of (ii) a 
hadron bubble in a bulk plasma. In both cases we normalize the free energy such that 
it is zero when there is no droplet resp. bubble, i.e. we calculate the free energy relative 
to an infinite hadron resp. plasma phase. In the adopted model, we have the following 
contributions to the free energy: Quarks (u, d, and s, and their anti-quarks), gluons, the 
bag contribution BVqgp, where Vqgp is the plasma volume, and the contribution from the 
three pions (using (|20| ) with p T (k, V) = 3^- and m = m w , this is 

F.(T,K) = ^ r^fcMna-e-v 7 ^/^), (52) 

where V n is the pion volume). The bag and pion contributions are thus simple and universal, 
the interesting part of the free energy is the quark and gluon contributions. 
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A. Free energy of a plasma droplet in a bulk hadronic medium 



The plasma droplet is an MIT bag with pions around it, so the calculation of the free 
energy is straightforward, i.e. we can calculate the energy levels of quarks and gluons in the 
bag, and perform the partition sum ([L9|) directly. This, of course, gives the true free energy. 

We can also use the MRE (p2|) to calculate the free energy (|20|). This way, two approx- 
imations are involved: (1) The spectrum is discrete, but the MRE treats the energy levels 
as continuous; and (2) we have discarded terms of the form 1/-R, 1/R 2 etc. in (|22|). Since 
we are considering the high temperature case, the first approximation is well justified: The 
low (E < T) energy levels, where the level spacing is large, do not contribute significantly 
to the free energy, whereas the main contribution to the true free energy (fL9D comes from 
the higher levels where the spacing is small. Thus, any difference in the free energy between 
a direct calculation and the MRE approximation is a measure of the importance of the 
neglected terms in fl2"2"| ) and/or the choice of truncation discussed in Sec. [TV]. 



1. Direct calculation 

In this case, we need to solve the set of equations an d then perform the sum (pj|) 

over these levels. We thus obtain the following equations (to be solved numerically) for the 
quarks (I — j ± 1/2) 

k 



j l (kR)=±— ji±i(kR), E = Jk* + m*. (53) 

E + m q v y 

For the TM-gluons 

3i(kR)=0, (54) 

and for the TE-gluons 

Jl (kR)(l + l) = kR Jl+1 (kR). (55) 

These equations provide a series of solutions, that we label Ej tn for the quarks and kf n , a = 
TM, TE for the gluons. The contribution to the free energy from quarks and gluons are 
then 

oo 

F q (T,R) = -QT £ (2j + l) £ JXl + e-WT) (5g) 

j=l/2,3/2,.. l=j±l/2 n 

and 

oo 

F g (T,R)=8Tj2(2l + l) E EMl-e^- /T ) (57) 

1=1 a=TM,TE n 

respectively. Note that when dealing with each angular momentum component separately, 
the degeneracy factors gi only account for the degeneracy due to color and particle/anti- 
particle, so we have g q = 6 for each quark flavor and g g = 8 for the gluons. 



12 



We imagine a spherical plasma droplet of volume Vqgp = ^-R 3 embedded in a large 
volume of pions, the total volume of this system being V^. The pions therefore inhabit a 
volume V n = Voo — Vqgp- Since we calculate the free energy relative to a system with no 
plasma droplet, and pions in the whole volume V^, the effective pion volume in Eq. (|52|) is 
—Vqgp- Summing all the contributions we obtain Fig. |j. In Fig. |5| we show the different 
contributions to Fig. |j. We shall comment on these figures later. 

2. Using the MRE 

Now let us see how the MRE approximation handles the plasma droplet. The difference 
from the sum over states method lies entirely in the calculation of the quark and gluon 
contributions. The pion and bag contributions are the same as before. Now, we use ( p0|) 
for the quarks and the gluons with the MRE density of states (p2|) . The result is shown in 
Fig. H Figure |7] shows the different contributions for T = 152 MeV. 

3. Comparison 

Comparing Fig. [| and Fig. |6| we see that although they agree qualitatively (except for 
R — ► 0, where the MRE is dominated by a curvature term proportional to R, whereas 
the sum over states for massless particles behaves as R 3 ), there are significant quantitative 
differences even at large radii. At first sight this is surprising. The MRE should be a good 
approximation, since each of the terms in ( f22"l) is derived analytically (albeit in the limit 
kR ^> 1). But the MRE is an approximation for the density of states, not for the free 
energy itself. Because the free energy is an integral over the density of states, it "picks up" 
the wrong behavior of the MRE at low energies, and "remembers" this error even at larger 
radii. This is why the free energy calculated using the MRE is not in quantitative agreement 
with the correct free energy in Fig. |], although the MRE is a good approximation for the 
density of states in the limit kR ^> 1. 

Looking at Fig. 0, we see that the main difference between the sum over states approach 
and the MRE approximation is due to the gluons. The gluon free energy being positive for 
small (R < 1.5 fm) radii, corresponds to the density of states being negative. This can also 
be seen directly from (^) and (|2"(3|). Hence the gluon density of states in the MRE is not only 
wrong, it is unphysical at small radii. This behavior is due to the way the MRE handles 
the surface corrections, namely through (P3|)-(P^). Taking more terms (oc l/R, oc 1/R 2 , 
etc.) into account in (^2|) might cure this. We propose another resolution of the problem in 
Sec. 0. 

A few words about the physics implied by Figs. [| and ^. There are two minima, at 
R = and at R = oo, separated by an energy barrier. The picture is therefore the following: 
When the temperature is T < To, no stable droplets can form. Even when T > Tq there 
is an energy barrier to pass before stable droplets of quark-gluon plasma can exist within 
the pion phase. But once a droplet is created (from fluctuations) with radius R > 2 fm, 
it will expand unimpeded. The energy barrier, of course, is due to the surface terms. The 
most important quantitative difference between Fig. [| and Fig. |6| is the height of the energy 
barrier separating the two minima at R = and R = oo. The height of the barrier is related 



13 



to the nucleation rate, in this case the plasma formation rate when heating a hadron gas 
e.g. in a heavy ion collision. Thus, although generally small, the surface contribution to the 
free energy has important implications. 



B. Free energy of a hadron bubble in bulk plasma 

This configuration is a vacuum bubble with hadrons (i.e. pions, in our model) inside, 
surrounded by plasma. In order to emphasize the essential points, we start by focusing on 
the surface free energy of quarks and gluons, since this is the interesting, non-universal, part 
of the free energy. 



1. Surface free energy of a vacuum bubble 

We shall compare results for this surface free energy as calculated by the MRE(— R) con- 
jecture and by the concentric spheres method. We shall also compare with results obtained 
by the phase shift approach. 

The figures [8], || and |10| show the surface free energy of massless quarks, massive quarks 
and gluons respectively. We show the results for just one temperature, T = 152 MeV, 
but the picture is qualitatively the same for other temperatures (T = 50,100,160 MeV). 
These figures show that the MRE(— R) works quite well in the case of quarks (massless and 
massive), but less well for gluons. The slope of the gluon MRE(— R) curve is correct, but 
there is an offset of about 1000 MeV (cf. Fig. 0). This is because the MEE(--R) for the 
density of states, although a good approximation at large values of kR, is wrong at small kR, 
and since the free energy at a given large radius is an integral over the density of states also 
at small momenta, the bad behavior of the MRE(— R) at small kR affects the free energy 
even at large radii (cf. remarks about the MRE for the plasma droplet). The same remarks 
apply to the quarks, but here the effect is much less pronounced. To support this picture 
we show in Fig. [11] a comparison of the quantity j^p(kR), p being the density of states of 
the eight gluons, calculated by the MRE(— R): j^Pmre(—R) {kR) = -^(kR) 2 + y, and by 



the phase shift approach (|33|). In contrast to other figures in this section, Fig. Ill] includes 



the volume contribution. Figure [12| compares all three methods of calculating the surface 
free energy, here shown in the case of gluons. We expect the concentric spheres method to 
be a correct way of calculating the free energy, as long as we stay in the regime i?Ci?2. 
The agreement of the concentric spheres method and the phase shift approach at R < 6 fm 
in Fig. [12| suggests that both approaches are valid ways to calculate the surface free energy 
and/or the density of states of quarks and gluons outside a vacuum bubble, for the phase 



shift approach presumably at any radius. Referring back to Fig. |TT] it is therefore clear that 
the MRE(— R) approach is inadequate at small kR. This was to be expected. First, the 
original MRE (for positive radii) is derived in the limit of large kR, and second, we have seen 
that the MRE is unphysical at small radii, so it is not surprising that also the MRE(— R) 
has problems in this regime. 

We thus conclude, that as an approximation of the density of states, the MRE(— R) for 
gluons works well at large values of kR, but is incorrect at small kR. We have only shown 
the gluon data, but a similar conclusion is valid for the quarks, although the error at small 
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radii is less important than in the case of gluons. Further, we have argued that as far as 
the free energy is concerned, we should be even more careful when applying the MRE(— R), 
as the bad behaviour of the MRE(— R) at small radii manifests itself as an error in the free 
energy even at large radii. Finally, our results obtained with the concentric spheres method 
are consistent with the phase shift formula (p7|), which seems to be an accurate way to 
calculate the density of states of quark-gluon plasma outside a vacuum bubble. 



2. The total free energy 



Knowing the contribution from the surface to the free energy of quarks and gluons outside 
the vacuum bubble, we can easily calculate the total free energy of the whole configuration: 
The volume contribution of the plasma is calculated using fl20|) with the smoothed density of 
states (^Tj) inserted. The pion contribution is given by fl52|), and the bag constant contributes 
a term BVqgp as always. Adding these contributions we obtain the results in Fig. [T3 . 



The interesting part of this figure is the minimum of the free energy at R = 1-2 fm for 
temperatures well above the transition temperature. Mardor and Svetitsky JTJ found a 
similar minimum in the free energy using the same model as described in this paper, but 
calculating the free energy of the plasma using the phase shift approach, whereas here we 
have applied the concentric spheres method. 



IV. CORRECTIONS TO THE MRE 



We have seen that the MRE and the MRE(— R) for gluons as it stands in (^), ( [25] ) and 
(p6|) have problems at small values of kR, leading to errors in the free energy even at large 
radii. Although we have numerical methods to calculate correctly the free energy both in 
the plasma droplet case and in the vacuum bubble case, we would like to be able to use some 
MRE approximation to gain physical insight, and for practical computations because the 
direct methods are numerically demanding. In this section, we investigate how to modify 
the MRE and the MRE(— R), in order for these approximations to describe more correctly 
the density of states of the configuration in question. 



A. Gluon droplet 



Previously in this paper, we have shown how to calculate exactly the free energy of 
abelian gluons in an MIT bag. That was, however, a numerical computation. Analytic 
calculations of the free energy, not using the MRE, have also appeared in the literature. 
Using the same model for the gluons as we do, De Francia |T6| , |T7|| finds for the difference 
AF = F„, uons — -fgiuons,MRE m the limit of large RT 

AF 5 

~Y = -0.874- ~ln(itT) + ... (58) 

where the dots indicate terms of higher order in (RT)' 1 . (De Francia gives such terms 
explicitly, but they are too small to be relevant in our analysis.) Note that (|58D is calculated 
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for one abelian gauge field, and should thus be multiplied by 8 in order to describe the gluon 
free energy. 

The main problem with the MRE is that it predicts a negative density of states at small 
kR where in reality there are no states, see Fig. [L4[ An error in the density of states at 
small values of kR is particularly severe, since here the statistical factor in the integrand of 
the free energy is large. In the following, we will show that using a reduced density of states 
of the form 

n (ub\ J°> 0<kR<A 

^PA(kR) = I 2 _ {kR)2 _ | ; kR> A ^ 

cures most of the problems of the MRE. We shall refer to this density of states, p\, as the 
MMRE (modified MRE), since it consists of the usual MRE-contributions for kR > A, but 
is truncated below kR = A. When we are in a regime where RT 3> 1, we can find an 
approximate analytical expression for the correction AF = F g i uonSj MMRE — -Pgiuons,MRE to the 
free energy induced by using p\ instead of p as the density of states of gluons: 

¥ = -16 (ln(A)(A A 3 _ _L A ) _ J_ A 3 + J_A + m ±A - A A 3)) . (60 ) 

We can fix the value of A by matching the coefficient of ln(i?T) to the analytical result of 
De Francia, i.e. solving 16(^A 3 — -J^A) = —5, which has A = 0.832 as the relevant solution, 
and the free energy (|60|) becomes 

AF 

— = -6.352 - 5 \n(RT). (61) 
The fact that 6.352 ~ 8 • 0.874 shows the consistency of the procedure, cf. Eq. (]58|) . 



In Fig. 15 we compare the proposal (B9f) with a direct calculation (summing over energy 



levels) of the free energy. Also shown is the free energy calculated using the usual MRE. 



B. Vacuum bubble 

Balian and Duplantier [18[] have calculated the Casimir energy of a perfectly conducting 
spherical shell. They find (in the large RT limit, and again not quoting terms of higher 
order in (RT)^ 1 ) 

AF = _0 : 769 _ ln(FT) 

T 4 4 [ ' 

In our language, 8 • AF is the sum of (i) the surface free energy of gluons inside an MIT 
bag, and (ii) the surface free energy of gluons outside a vacuum bubble. Using this and De 
Francia's calculation ([58"D, we can deduce the corrections to the MRE(— R). We obtain for 
the difference AF vac = F g i uonSiCorrect ed — F g i uonSi MRE(-R) 

AF vac = T(5.454 + 31n(FT)). (63) 

Like in the gluon droplet case, we can advise a modification of the MRE(— R), which works 
very well. Specifically, we propose the following: The reduced density of states of gluons 
outside a vacuum bubble, is 
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ix ,, m _jO, 0<kR< 0.458 , . 

rP™W K ) ~ \ -l{kRf + 1, kR> 0.458 lb4j 

where > is the radius of the vacuum bubble. We refer to this density of states as 
the MMRE(— i?). The value of kR = 0.458 where we cut the MRE(--R) is fixed by the 
same procedure as in the gluon droplet case. The difference AF = -F g i uonSj MMRE(-R) — 

-Pgluons,MRE(-R) I s then 

AF = T(5.417 + 3 \n(RT)), (65) 

showing the consistency of the procedure (cf. Eq. (|63|)). The fact that we should cut the 
density of states at kR = 0.458 and not kR = 0.832 as in the gluon droplet case, reflects 
the asymmetry between the gluon droplet and the vacuum bubble configurations. 

In Fig. |Tj] we compare the different methods of calculating the free energy of gluons 
outside a vacuum bubble. The simple MMRE(— R) suggestion is in nice agreement with the 
phase shift approach, which (based on our calculations in the previous sections) we consider 
the most accurate way of calculating the free energy. 



V. CONCLUSION 

This paper has a twofold purpose. First, we have introduced the concentric spheres 
method as a way to calculate the free energy of quark-gluon plasma outside a pion bub- 
ble, confirming the peculiar results of Mardor and Svetitsky PJ, that, within the MIT bag 
model, this free energy has a minimum at non-zero radius even well above the transition 
temperature. 

Second, we have shown that terms beyond volume, surface, and curvature are necessary 
in order to reproduce the free energy of plasma droplets and vacuum bubbles within the 
multiple reflection expansion, especially for the gluon contributions. We have discussed the 
reasons for this, and based on previous calculations lfL6l,|18| , we extract correction terms to 
the free energy, which can be understood from a physically motivated truncation of the 
density of states. 

Our calculations are all performed in the limit of vanishing chemical potentials. The 
results are thus relevant to investigations of the cosmological quark-hadron transition, and 
possibly to forthcoming ultrarelativistic heavy ion collision experiments at RHIC and LHC. 
While these are certainly interesting prospects, we plan to extend our analysis to situations 
of finite chemical potential, such calculations being relevant to a wider range of applications 
including e.g. neutron stars. 
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FIGURES 




FIG. 1. Left: A plasma droplet with non-perturbative vacuum outside; this is essentially 
an MIT bag. Right: A vacuum bubble surrounded by plasma, boundary conditions being those 
of the MIT bag, but corresponding to no flux of plasma into the bubble. The phase outside a 
droplet/bubble extends to infinity. 




FIG. 2. The concentric spheres configuration: The inner sphere has radius R±, the outer 
sphere has radius R 2 - There is non-perturbative vacuum inside the inner sphere and outside 
the outer sphere. Quark-gluon plasma is confined between the two spheres by MIT bag bound- 
ary conditions corresponding to no flux of plasma across the spheres. At the outer sphere the 
boundary conditions are the usual MIT bag conditions ([|)-([]), but at the inner sphere we use 
nV( x )\ r=Rl = -n»{x)\ T=R2 = (0,f/|x|). 




FIG. 3. How to extract the surface free energy in the concentric spheres method: Shaded areas 
are plasma, white areas are non-perturbative vacuum. Contributions to free energy (left to right): 
F tot ai(Ri,R2), F tota i(R 2 ), F vo i ume {Ri) and F surface (Ri). F tot ai(Ri, R2) and F tota i{R2) are calcu- 
lated by summation over energy levels of the particles in the relevant configuration. F vo i ume (R\) 
is only the volume free energy of particles occupying a volume ^-Rf. The purpose is to calculate 
the contribution to the free energy from the inner surface, F sur f ace (Ri), and this can be done as 
in @. 
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1 2 3 4 5 
R(fm) 

FIG. 4. Total free energy (calculated directly by summation over states) of a quark-gluon 
plasma droplet of radius R surrounded by pions. The phase transition temperature is set to 
Tq = 150 MeV. Results are shown for several temperatures around Tq. 
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FIG. 6. As Fig. ||, except that the free energies are calculated using the MRE approximation. 
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FIG. 8. A comparison between two different methods of calculating the surface free energy 
of massless quarks outside a sphere (a "vacuum bubble") of radius R: (1) Multiple reflection 
expansion with the sign of R reversed (MRE(— R)), and (2) the concentric spheres method with 
an outer radius R2 = 20 fm. When R < R%/2 the two methods yield similar results. This 
suggests two things: (A) When R <C R2 the interactions at the outer surface are unimportant, and 
(B) MRE(— R) describes adequately the way the inner surface alters the density of states in the 
case of massless quarks. The temperature is T = 152 MeV. 
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FIG. 9. Same as Fig. ||, but using a quark mass of 150 MeV. In this case, the effect of the 
outer surface is not visible until R is quite close to the outer surface at R2 = 20 fm. Again, the 
MRE(— R) seems to be a satisfactory description. 
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FIG. 10. Same as Fig. ^ but for gluons. The MRE(— R) is not as good a description as in the 
case of the quarks. 
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FIG. 11. The "reduced density of states" ^p(kR) of gluons outside a vacuum bubble, calcu- 
lated by the phase shift approach, and using the MRE(— R). The phase shift method yields the 
more correct result (cf. Fig. |l2|). Note that this is the density of states relative to a situation with 
no vacuum bubble, so p < in this case is not unphysical. 
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FIG. 12. Surface free energy of gluons outside a vacuum bubble of radius R calculated in the 
three different ways described in the text. The outer radius used in the concentric spheres method 
is i?2 = 20 fm. The phase shift approach agrees with the concentric spheres method for R < 6 fm, 
whereas the deviation between these two methods at R > 6 fm is due to the influence of the outer 
surface on the result of concentric spheres method. To emphasize the differences between the three 
methods, we have subtracted the volume free energy from the phase shift results, so that only the 
surface contributions are shown in this figure. The temperature is T = 152 MeV. 
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FIG. 13. Free energy of a pion bubble surrounded by quark-gluon plasma, normalized so that 
the free energy of a pure plasma without pions is zero. Curves are shown for temperatures above 
and below the transition temperature, To = 150 MeV. The surface contributions from the quarks 
and gluons are calculated by the concentric spheres method with an outer surface of R2 = 20 fm. 
The minimum at R ~ 1 — 2 fm shows that in this model, bubbles of pions of this radius will form 
even for temperatures above To. Similar results were obtained by Mardor and Svetitsky using 
the phase shift method. 
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FIG. 14. Number of gluon states with energy less than k in an MIT bag of radius R, calculated 
(i) directly, solving ( |54| ) and (|55|) (discontinuous line, true values) and (ii) by the MRE (continuous 
line, approximation). Note that the MRE predicts a negative density of states at small values of 
kR. 
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FIG. 15. Gluon contribution to the free energy (in units of the temperature, T) of a plasma 
droplet of radius R, as a function of the dimensionless parameter RT. The unphysical behavior 
of the usual MRE ((p2|), (B^) and (^6|)) at small radii causes the free energy to deviate from the 
true free energy (calculated by summing over energy levels) even at large radii. Using the modified 
MRE, po.832> in ( p0| ) makes the free energy agree remarkably well with the true free energy. 
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FIG. 16. The free energy of gluons outside a vacuum bubble of radius R (normalized to the 
temperature, T) calculated in three different ways: (1) By the phase shift approach, which we 
consider an accurate procedure, (2) by the MRE(— R), and (3) using the MMRE(— R). 
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